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Abstract  K 

Ai r/sea  measurements  from  the  Long-^errn  Upper  Ocean  ^tudy  (LOTUS)  buoy 
in  the  Sargasso  Sea  arej[^analyzed  to  learn  how  the  diurnal  response  of  sea 

surface  temperature,  iaT\,  is  related  to  the  surface  heating,  H,  and  the 

■  ' 

wind  stress,  S.  Data  are  taken  from  the' LOTUS-3  and  LOTUS-5  records  which 
span  the  summers  of  1982  and  1983.  The  b^ic  data  are  shown  in  monthly 
plots,  and  tne  analyzed  daily  values  of  [aT^,  H,  and  S  are  given  in 
tables  and  in  figures. 

Analyzed  data  show  a  clear  trend  of  aT\  increasing  with  H  and 
decreasing  with  S.  A  best-fit,  three-parameter,  empirical  function  can  ac¬ 
count  for  90/'percent'of  the  variance  in  a  screened  subset  of  the  LOTUS  data 
(172  days)  and  81 ‘percent’ of  the  variance  of  tne  full  data  set  (361  days). 

The  analyzed  data  are  also  compared  with  a  theoretical  model  function 
now  used  for  ocean  predictions  in  the  Diurnal  Ocean  Surface  Layer  model 
(UOSL)  of  Fleet  Numerical  Oceanography  Center.  The  DOSL  model  function  was 
derived  from  the  assumption  that  wind-mixing  occurs  by  a  mechanism  of  shear 
flow  instability.  It  is  fully  predictive  and  shows  a  parameter  dependence 
consistent  with  the  LOTUS  data  over  a  wide  range  of  H  and  S.  The  DOSL 
model  function  can  account  for  almost  as  much  variance  as  the  best-fit  empir¬ 
ical  function.  F  m  A',  .'m.'  ^  .  iMr-ir'.'r  'i  ,  ; 


TABLE  OF  CONTENTS 


Page 


Abstract  .  i 

List  of  Tables  .  4 

List  of  Figures .  4 

1.  Introduction  .  b 

II.  LOTUS  Data  and  Analysis  Methods  .  b 

A.  Sea  Surface  Temperature  Response  .  b 

B.  Dai ly  Averages  .  6 

C.  Insolation  .  7 

1)  Maximum  Observed  Insolation  .  7 

2)  Integrated  Insolation  .  7 

3)  Least  Squares  Fit  to  Insolation  .  8 

0.  Horizontal  Advection  .  8 

III.  Results  9 

IV.  Fit  to  and  Comparison  with  Models .  14 

A.  Empirical  Model .  14 

B.  Positive  Bias  of  Analyzed  aT^  . lb 

C.  Theoretical  Model  .  lb 

V.  Conclusions  .  20 

Acknowledgments  .  22 

References  23 

Appendix  A,  Monthly  Time-Series  Plots  of  LOTUS  Data  .  24 

Appendix  B,  Listing  of  Analyzed  Data  .  38 

Appendix  C,  FORTRAN  Listing  of  DOSL  Model  Function  .  48 


LIST  UP  TABLES 


I.  Statistics  of  Data  /  Model  Comparisons 


LIST  OF  FIGURES 

1.  Analyzed,  full  data  set  . 

data  set  with  advection  periods  omitted  . 

3.  Screened  data  set  . 

4.  Screened  data  and  best-fit  model  function  . 

5.  Sea  surface  temperature  response 

at  very  low  values  of  heating  . 

6.  Screened  data  set  and  DOSL  model  function  . 

7.  Residuals  from  the  screened  data  set  and  OUSL  model  function 
A.l- 

A.ii  Monthly  time-series  plots  of  LOTUS  data  . 


I.  INTRODUCTION 

In  this  study  we  use  oceanographic  and  meteorological  field  observa¬ 
tions  to  examine  the  response  of  sea  surface  temperature  to  diurnal  heating 
and  wind  stress.  The  data  were  acquired  by  the  ]^ng-Term  IJpper  Ocean  ^tudy 
(LOTUS)  buoy  which  was  deployed  in  the  northwestern  Sargasso  Sea  (34’N,  70“W) 
from  198^!  to  1984  (Briscoe  and  Weller,  1984;  Deser,  Weller  and  Briscoe,  1983; 
Tarbell,  Montgomery  and  Briscoe,  1985).  The  LOTUS  data  are  well  suited  for 
this  study  because  the  signal  of  local  heat  storage  is  generally  large  com¬ 
pared  to  the  effects  of  horizontal  advection,  and  because  the  LOTUS  data 
provide  high  accuracy  and  temporal  resolution. 

The  aim  of  the  analysis  is,  first,  to  reduce  the  time  series  data  to 
daily  values  of  the  heating,  H,  wind  stress,  S,  and  surface  temperature 
response,  aT^,  and  second,  to  determine  the  dependence  of  aT^  upon  H 
and  S.  There  is  a  practical  value  to  this  result  —  the  function  aT^(H,S) 
defined  or  verified  by  these  data  can  be  used  to  forecast  or  hindcast  the 
diurnal  cycle,  and  a  scientific  value  —  these  data  show  how  a  rotating  fluid 
responds  to  a  stabilizing  heat  flux  and  an  imposed  stress. 

The  primary  purpose  of  this  technical  report  is  to  provide  time  series 
plots  of  the  original  LOTUS  data,  a  tabulation  of  the  analyzed  data,  and  a 
listing  of  the  FORTRAN  program  which  evaluates  a  theoretical  model  function. 

II.  LOTUS  OATA  AND  ANALYSIS  METHODS 

The  LOTUS  data  used  for  this  analysis  are  time  series  measurements  of 
wind  speed,  insolation  and  ocean  temperatures  at  0.6,  5,  10,  15,  50  and 

7b  m  at  15  min  intervals  (Deser  et^.,  1983;  Tarbell,  Pennington  and  Bris¬ 
coe,  1984).  Estimates  of  the  wind  stress,  r,  and  surface  heat  flux,  Q  = 

I  *  L,  where  I  is  insolation  and  L  is  the  heat  loss,  were  made  as  de¬ 
scribed  by  Stramma  ^  (l!J86).  The  data  sets  analyzed  are  from  the  de¬ 

ployments  LOTUS-3  (May  14,  1982  through  October  20,  1982)  and  LOTUS-5  (April 
15,  1983  through  October  30,  1983).  Monthly  plots  of  the  full  data  set  are 
in  Appendix  A.  In  Appendix  B  we  list  the  analyzed  daily  values  of  the  wind 
stress,  wind  speed,  insolation,  etc.,  so  that  other  investigators  will  have 
ready  access  to  our  intermediate  results. 


A.  Sea  Surface  Temperature  Response 

The  ocean's  response  to  the  forcing  of  the  winds  and  solar  heating  is 
characterized  here  by  the  amplitude  of  the  sea  surface  temperature  response, 
aT^,  where  T^  is  the  LOTUS  temperature  at  0.6  m  depth,  and  a  indicates 
the  change  (increase)  over  a  day.  If  tne  diurnal  cycle  is  the  dominant  sig¬ 
nal  in  T  ,  then  the  lowest  T  of  the  day  will  occur  between  0500  and 
s  s 

0900  hours  (all  times  quoted  in  the  text  are  local  solar;  local  noon  is  1730 
Z),  and  the  highest  T^  will  occur  between  1100  and  1800.  We  have  esti¬ 
mated  the  diurnal  response  of  sea  surface  by  subtracting  the  lowest  tempera¬ 
ture  Observed  between  0500  and  0900  from  the  highest  temperature  observed 
between  1100  and  1800.  The  result  is  listed  in  Appendix  B  as  "DEL  T". 

This  process  of  subtracting  the  lowest  from  the  highest  temperature 
will  alias  high  frequency  noise  (compared  to  diurnal)  into  larger  aT^.  In 
Section  I\/.B  we  will  evaluate  the  resulting  bias  in  tne  analyzed  aT^  by 
inspection  of  the  estimates  available  at  very  low  values  of  surface  heating. 

6.  Oai ly  Averages 

In  tne  following  analysis  we  compute  a  "daily"  value  as  an  average  over 
the  period  during  each  day  when  Q  >  0.  Let  tj^  be  the  time  when  Q 
first  becomes  positive,  and  let  t^  be  the  time  when  Q  becomes  negative 
at  sunset.  Daily  average  is  defined  here  as 


)  =  jpr- 

^  Q 


)  dt  , 


where  =  1/2  (t^  -  t^ ) .  (On  some  days  Q  may  become  negative  within  the 
time  interval  t^  <  t  <  t2  when  increased  cloud  cover  causes  the  insolation 
to  drop  below  the  value  of  the  heat  loss.  A  second,  alternate,  time  scale, 
PQ^,  (not  used  here)  was  defined  to  be  literally  the  half  time  that  Q 
was  positive  for  each  day.) 

Daily  average  values  are  tabulated  in  Appendix  B  as: 

"U",  daily  average  wind  speed  (m  s”^), 

"TAU",  daily  average  wind  stress  (Pa), 

_2 

"L",  average  heat  loss  (W  m  ). 
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standard  deviations  are  also  listed  as  “SD  U",  "SD  TALI',  and  "SD  L".  In  com¬ 
puting  SD  TAD  we  sum  the  standard  deviations  of  the  east  and  north  components 
to  take  account  of  varying  wind  direction.  In  some  of  the  later  comparisons 
with  models  we  will  screen  out  those  days  having  highly  variable  or  irregular 
forcing  and  retain  the  days  which  have  comparatively  steady  forcing. 


C.  Insolation 

The  amplitude  of  insolation  has  been  estimated  in  three  ways. 

1)  Maximum  Observed  Insolation 

The  maximum  insolation  observed  on  each  day  was  extracted  from  the  rec¬ 
ords,  and  is  listed  as  I^.  If  there  were  no  clouds,  then  1^^  together  with 
the  duration  of  the  insolation,  D,  would  completely  characterize  the  inso¬ 
lation  for  any  day.  However,  cloud-free  days  are  the  exception  rather  than 
the  rule,  and  we  have  therefore  defined  two  other  more  useful  measures  of 
insolation. 

2)  Integrated  Insolation 

The  integral  of  the  insolation  has  been  used  to  define  I2  as 


f  M 

t. 


I  dt 


'2  = 


Dr 


where 


I  =  observed  insolation. 


t^  =  the  start  of  insolation. 


t^  =  the  end  of  insolation. 


D  =  duration  of  the  insolation, 

rt/i 


^4  "  ^3’ 


r  = 


1 


OF. 


m 


Jt 


F(t,  i),  o)dt  =  0.57,  a  constant. 


F(t,  0,  a)  =  theoretical  clear  sky  insolation;  function  of  year  day,  t, 
latitude  i  =  34’n,  and  clear  sky  transmittance,  a  =  0.8. 


F^  is  the  daily  maximum  (List,  1958). 


Because  the  insolation  records  have  some  slight  noise,  t,  was  esti- 

-2  ^ 

mated  to  occur  a  half  hour  before  I  >  30  W  m  ,  and  similarly  t.  was  ta- 

-2  ^ 

ken  to  be  a  half  hour  beyond  I  £  30  W  m  in  late  afternoon.  The  error  in 
defining  the  period  of  insolation  in  this  manner  is  ±  15  minutes  for  most  days, 
which  gives  about  4  percent  error  in  the  computed  amplitude  (I^)  of  the 
insolation. 

3)  Least  Squares  Fit  to  Insolation 

As  an  alternative  measure  of  insolation,  we  have  also  carried  out  a 
least  squares  fit  of  the  normalized  theoretical,  insolation  function  F  onto 
observed  insolation  to  compute  the  amplitude  (Appendix  B).  The  stan¬ 
dard  deviation  of  the  fit,  listed  as  "SD  provides  a  convenient,  objec¬ 

tive  measure  of  the  variability  of  insolation  due  to  cloud  cover.  In  a  later 
analysis  we  screen  out  days  with  highly  variable  insolation  by  setting  an 
upper  limit  on  the  ratio,  SD  13/12* 

D.  Horizontal  Advection 

The  analysis  thus  far  has  been  concerned  only  with  the  vertical  fluxes 
of  heat  and  momentum.  Our  implicit  assumption  of  a  local  balance  is  valid 
as  long  as  there  is  no  horizontal  advection  occurring  in  the  water  column. 
However,  the  effects  of  horizontal  advection  are  evident  in  the  LOTUS  data 
as  occasions  when  the  water  column  (above  75  m)  underwent  depth-independen.t 
temperature  changes.  Since  we  were  unable  to  account  for  the  effects  of 
horizontal  advection  explicitly,  we  sought  to  identify  at  least  the  most 
obvious  occurrences  of  horizontal  advection  so  that  we  could  eliminate  the 
corresponding  days  from  the  analysis. 

A  local  heat  balance- analysis  was  performed  by  comparing  the  net  sur¬ 
face  heat  flux  with  the  observed  heat  storage  in  the  water  column. 


where  B  is  the  observed  heat  storage  to  25  m  depth,  and  aB  =  B(t2)  - 
B{c,).  The  values  of  EPS  are  listed  in  Appendix  B.  This  alone  does 
not  effectively  identify  horizontal  advection  since  vertical  advection  can 
also  cause  significant  changes  in  heat  storage  without  changing  T^  (note 
especially  the  period  10  to  30  September  1983,  Figure  A. 12,  when  there  were 
very  large  changes  in  temperature  at  25  m  and  deeper,  but  no  corresponding 
changes  in  the  near  surface).  We  have  therefore  used  the  tabulated  EPS 
values  as  a  guide,  and  proceeded  to  subjectively  identify  periods  where 
horizontal  advection  seemed  to  be  important.  These  were  noted  by  setting 
the  flag  A  =  1  (Appendix  B),  and  included  only  44  days  (two  periods  in 
1982  and  three  in  1983),  or  about  12  percent  of  the  complete  data  set. 


III.  RESULTS 

Our  analysis  of  the  full  data  set  is  presented  graphically  in  Figure  1 
where  the  diurnal  response  of  sea  surface  temperature  is  plotted  against 
independent  variables  proportional  to  the  heating  and  wind  stress.  The 
"heating"  variable  was  computed  as: 

Q  Pq 

H  = - ^  (C  m) 

D  U 
0  p 

where 

Q  =  +  L  (W  m“^), 

=  1023  (kg  m~^) ,  and 
Cp  =  4183  (J  kg"^  C"^). 

The  product  Q  Pq  is  proportional  to  the  (warming)  heat  flux  supplied  to  tne 
ocean;  division  oy  Cp  gives  kinematic  units,  C  m,  which  are  more 
readily  interpreted.  In  a  similar  way  the  "stress"  variaole  was  computed  as: 

TP 

b  =  - ^  (m^  s"^), 

"o 

where 

I  =  daily  average  wind  stress  (Pa), 

1/2 

P  =  acceleration  time  scale  =  l/f[2  -  2  cos(fPf,)]  (s),  and 

'  -1  ^ 
f  =  Coriolis  parameter  (s  ). 


I], 


In  this  final  regard  the  scaling  begins  to  follow  the  theory  developed 
by  Price,  Weller  and  Pinkel  (1986;  hereafter,  PWP)  where  the  vertical  mixing 
due  to  wind  stress  was  presumed  to  occur  through  shear  flow  instability 
caused  by  wind-driven  currents.  [The  amplitude  of  the  wind-driven  current 
(or  diurnal  jet)  was  estimated  to  be  W/D  where  D  is  the  trapping  depth, 
and  ~  H/D.]  The  f-dependence  of  P^  takes  account  of  rotation  which 
in  this  case  is  small  enough  that  P^  =  P^.  In  effect,  the  heating  rate 
Q  and  the  wind  stress  ^  are  multiplied  oy  nearly  the  same  time  scale 
(which  does  vary  from  day  to  day). 

While  there  is  a  good  deal  of  scatter  in  the  data  of  Figure  1,  it  is 
also  apparent  that  there  is  a  significant  functional  dependence  of  aT^ 
upon  H  and  S.  Within  the  variable  range  sampled  here,  aT^  increases 
with  increasing  H  and  with  decreasing  S  —  the  qualitative  result 
expected. 

Before  we  attempt  to  define  or  test  a  model  function,  we  first  screen 
the  data  to  omit  the  days  which  one  would  expect,  a  priori ,  to  be  unsuitable 
for  defining  a  function  aT^(H,S).  That  is,  we  assume  that  if  the  wind 
was  steady  during  the  day  and  H  was  regular  (not  intermittent  due  to  cloud 
cover),  then  some  function  aT^(H,S)  should  obtain.  On  the  other  hand, 
if  the  wind,  say,  were  highly  variable  during  a  day,  then  there  are  addi¬ 
tional  degrees  of  freedom  present,  and  a  function  dependent  upon  S  alone 
might  not  be  appropriate.  The  LOTUS  data  set  is  large  enough  that  we  can 
omit  days  with  irregular  or  variable  forcing  and  still  retain  a  large  sample 
spanning  a  wide  parameter  range. 

In  Figure  2  we  have  omitted  the  days  identified  as  showing  some  effects 
of  horizontal  advection  (discussed  in  Section  II. D),  and  in  Figure  3  we  have 
also  omitted  days  having  either  highly  variable  insolation,  SO  I3/I2  > 

1/2,  or  highly  variable  wind  stress,  SO  TAU/TAU  1/2.  The  result  of  this 
screening  is  to  eliminate  some  of  the  points  which  lie  furthest  from  the  mean 
trend  of  the  data.  The  last  version.  Figure  3,  will  be  used  for  model  defi¬ 
nition  and  testing,  but  we  also  make  comparisons  with  the  full  data  set. 
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IV.  FIT  TO  AND  COMPARISON  WITH  MODELS 

These  data  may  be  used  to  define  purely  empirical  model  functions  and 
to  test  theoretical  model  functions.  To  quantify  data/model  comparisons,  we 
have  calculated  statistics  of  the  deviations. 


T'  =  aT^  -  aT^ 
s  m 


where  aT^  is  aT  evaluated  by  the  best  fit  or  theoretical  model  at  the 
H,  S  of  the  corresponding,  observed  aT^.  The  ensemble  Average  Deviation 

IS 

> 

AD  =  <T'>  =  ^  I  T‘  , 


where  N  is  the  number  of  days; 

the  Percentage  of  Variance  accounted  for  by  the  model  is 


PV  =  100 


1  - 


t.2 

<T'  > 


and  the  Correlation  Coefficient  is 


CC  = 


<aT„  aT^> 
m  s 

~~rTrr~7~rm 

<aT  >  <aT^^> 

m  s 


where  (  )  indicates  departure  from  the  ensemble  mean.  (PV  and  CC  differ 
in  that  CC  is  independent  of  bias  error.  Here  the  bias  error  is  small  com¬ 
pared  to  random  errors,  and  we  emphasize  PV  in  our  discussion.) 


A.  Empirical  Model 

A  purely  empirical  model  function  is  defined  from  the  data  by  maximiz¬ 
ing  PV  for  the  three  parameter  function, 

aT{H,S)  =  aH®  exp(-S/Y). 
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This  functional  form  was  chosen  by  inspection  of  the  data,  and  under  the  as¬ 
sumption  that  aT  should  vanish  as  H  >  0,  and  should  be  finite  as  S  >  0. 
Values  of  o,  6,  Y  which  maximize  PV  were  found  by  a  searching  method  to 
be: 

o  =  0.20  ±  0.03  (C)  , 

B  =  1.40  *  0.1  and 
Y  =  0.80  0.05  (Pa). 

These  give  PV  =  90  (see  also  Table  I),  showing  that  there  is  indeed  a  strong 
dependence  of  aT^  upon  the  presumed  independent  variables,  H  and  S. 

This  best  fit  function  is  plotted  as  a  surface  along  with  the  screened  data 
in  Figure  4. 


8.  Positive  Bias  of  aT„ 

s 

A  consistent  result  of  the  liata/model  comparisons  is  that  AD  tends 
to  be  positive,  i.e.,  on  average  the  observed  aT^  lie  slightly  above  the 
model  prediction  (as  can  be  seen  in  Figure  4).  This  is  at  least  partially  a 
result  of  a  bias  error  inherent  in  the  day-night  differencing  procedure  used 
to  estimate  the  diurnal  response  of  surface  temperature  (Section  II.A).  To 
checx  this  we  have  plotted  the  aT^  available  at  small  values  of  H, 

Figure  5.  Note  that  aT^  tends  to  remain  slightly  positive  as  H  goes 
to  zero;  <aT^>  =  0.07  C  for  H  <  0.05  C  m,  which  is  unphysical,  and 
presumed  to  be  a  bias  error  of  the  analysis. 

C.  Theoretical  Model 

A  theoretical  model  function  derived  by  PWP  and  now  in  use  at  Fleet 
Numerical  Oceanography  Center  as  part  of  the  Diurnal  Ocean  Surface  Layer 
(DOSL)  model  may  be  tested  using  these  data.  The  only  additional  parameters 
needed  to  evaluate  the  model  (see  Appendix  C  for  model  listing)  are  those 
that  define  the  optical  properties  of  the  water.  For  the  LOTUS  site  we  use 
Type  I  parameters  appropriate  for  very  clear  ocean  waters  (Stramma  et  al . , 
1986;  Paulson  and  Simpson,  1977:  long  wave  extinction  scale,  b^  =  0.35  m; 
short  wave  extinction  scale,  82  =  23  m,  and  the  fraction  of  long  wave 
insolation,  R  =  0.58) . 

The  statistics  of  this  model/data  comparison  are  in  Table  I,  and  the 
result  is  plotted  in  Figure  6.  Note  that  AD  is  again  slightly  positive  and 
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Table  I 


Statistics  of  Data  /  Model  Comparisons 


Best  Fit, 
Empirical  Model 


DOSL  Model 


Jv'Cv 

V’'  -  ' 

Data  Set 

Number  of 

Days 

361 

Average 

•-V>’ 

Deviation, 

0.06 

C 

Percent 

Variance 

81 

o 

Cross- 

m 

Correlation 

0.82 

Screened 
Data  Set 


Full 

Data  Set 


Screened 
Data  Set 


Screened  data  set  and  the  best-fit  empirical  model  function  shown 
as  a  surface.  Data  points  which  lie  above  the  surface  are  solid, 
those  below  are  open.  Vertical  lines  connect  the  points  to  the 
model  surface. 


HEATING,  C  m 


Sea  surface  temperature  response  at  very  low  values  of  heating. 
Note  that  the  estimates  have  a  small  mean  value,  ~  0.07  C,  as 
heating  vanishes,  probably  because  of  a  bias  error  inherent  in  the 
day/night  differencing  method  used  to  define  the  diurnal  response. 


that  PV  =  87,  or  nearly  as  successful  as  the  best-fit  empirical  function. 

The  empirical  model  function  and  the  OOSL  model  function  differ  signifi¬ 
cantly  only  in  the  limit  of  very  large  H  and  vanishing  S  where  there  are 
not  enough  data  to  tell  which  is  better. 

The  DOSL  model  function  appears  to  follow  the  trend  of  the  data  reason¬ 
ably  well  (suggested  also  by  large  PV).  To  check  explicitly  whether  there 

may  be  coherent  structure  in  the  deviations,  we  have  smoothed  and  interpo¬ 
lated  T'  to  a  regular  grid,  and  plotted  the  resulting  local  average  T' 

(local  in  H,S)  as  a  surface  in  Figure  7  left.  This  surface  lies  slightly 

above  0  almost  everywhere  (by  about  0.05  to  0.1  C),  consistent  with  the 
bias  noted  before.  There  is  no  obvious  low  mode  structure  to  the  surface, 
which  suggests  that  the  model  function  dependence  upon  H  and  S  is  reason¬ 
ably  consistent  with  the  H,  S  dependence  of  the  LOTUS  data.  Said  differ¬ 
ently,  the  deviations  do  indeed  appear  to  be  random  (aside  from  the  bias  of 
the  analysis). 

In  a  similar  way  the  local  average  root  mean  square  deviation  was  com¬ 
puted  and  plotted  in  Figure  7  right.  This  surface  does  have  a  significant 
low  mode  structure  with  the  rms  T'  increasing  as  does  aT^.  In  the  range 
where  aT^  =  2  C,  the  rms  T'  =  0.4  C.  (In  the  limit  of  small  H,  large 
S,  the  rms  T'  =  0.1  C,  which  is  essentially  the  statistics  of  the  data 
itself.)  A  random  error  which  increases  with  the  signal  might  result  from 
variable  wind  stress  or  insolation  (recall  that  SO  TAU/TAU  is  allowed  to  be 
as  large  as  1/2  here). 

A  statistical  comparison  of  the  DOSL  model  with  the  full  data  set  shows 
that  there  is  still  a  significant  dependence  of  aT^  upon  H  and  S  (PV  = 
74,  Table  I),  though  on  some  specific  days  the  deviation  may  be  quite  large. 

V.  CONCLUSIONS 

The  LOTUS-J  and  LOTUS-5  data  sets  have  been  used  to  define  and  test 
model  functions  that  relate  the  diurnal  response  of  sea  surface  temperature 
to  tne  imposed  surface  heat  flux  and  wind  stress.  Either  a  best  fit  empiri¬ 
cal  function  or  the  DOSL  model  function  developed  by  PWP  can  account  for 
roughly  80 »  of  the  aT^  variance  observed  in  the  LOTUS  data  set.  The  DOSL 
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Figure  7:  The  average  (top)  and  root  mean  square  (bottom)  field  of  the  devi 
ations,  T',  from  Figure  6.  The  surfaces  were  constructed  by 
smoothing  over  neighboring  values.  The  heavy  border  shows  where 
deviation  =  0. 

Note  that  the  average  value  tends  to  be  about  0.1  C 
throughout  the  full  range  of  H  and  S,  with  little  evidence  of 
low  mode  structure.  On  the  other  hand,  the  rms  value  is  largest 
where  aTj  is  also  largest.  This  suggests  that  there  may  be  a 
"hidden  variable"  not  accounted  for  by  the  present  scaling  and 
model  assumptions  (e.g.,  using  only  the  daily  mean  value  to  rep¬ 
resent  the  wind  stress). 


function  has  the  advantage  of  showing  explicit  dependence  upon  sea  water  op¬ 
tical  properties  and  latitude,  while  in  the  empirical  function  this  depen¬ 
dence  is  implicit  in  the  parameter  values. 

The  DOSL  model  function  was  based  upon  the  idea  that  the  currents  gen¬ 
erated  by  wind  stress  cause  vertical  mixing  by  a  mechanism  of  shear  flow  in¬ 
stability.  The  resulting  function  aT^IH.S)  which  follows  from  this  idea 
has  a  parameter  dependence  consistent  with  the  LOTUS  field  data.  This  is 
indirect  evidence  that  the  mixing  assumption  was  a  good  one  (but  to  make  a 
strong  test  requires  sensitivity  testing  and  alternate  model  testing  not  at¬ 
tempted  here). 

The  tests  of  the  OOSL  model  function  made  here  do  suggest  that  that 
function  is  appropriate  physically.  That  is,  if  given  accurate  forecasts  of 
heating  and  wind  stress,  the  DOSL  model  should  return  an  accurate  upper  ocean 
forecast.  Of  course,  on  any  specific  site  or  day,  the  ocean's  response  could 
be  dominated  more  by  advection,  or  by  sub-mesoscale  variability  of  the  winds, 
than  by  the  forecast  diurnal  response  (recall  the  larger  scatter  in  the  full 
data  set  compared  to  the  screened  data  set.)  It  may  be  that  the  large  scale 
(atmospheric  mesoscale)  patterns  of  diurnal  warming  will  be  forecast  more 
successfully  than  will  the  point-wise  response  studied  here. 

The  most  readily  observable  quantities  are  the  cloud  cover  (needed  to 
calculate  H)  and  the  sea  surface  temperature  itself.  This  suggests  that  an 
inversion  of  aT^(H,S)  to  estimate  S  (or  wind  speed)  might  be  useful  at 
least  in  the  low  S  regime  where  there  is  some  sensitivity.  As  a  qualita¬ 
tive  example,  Stramma  (1986)  have  shown  that  regions  of  large  aT^ 

are  coincident  with  regions  of  weak  sea  surface  pressure  gradient,  e.g., 
ridges  of  the  marine  high  pressure  systems.  Thus,  given  a  map  of  H  and  aT^ 
a  forecaster,  or  an  analysis  program,  could  easily  sketch  in  the  locations 
and  perhaps  some  measure  of  the  width  of  low  wind  speed  regions  over  the 
world  ocean. 
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F igure  A. L :  Monthly  time  series  from  May,  1982.  The  periods  during  which 
~  advection  appears  to  be  important  (discussed  in  Section  II. D) 

are  denoted  by  the  bar  labelled  A. 
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Figure  A .4:  Monthly  time  series  from  August,  1982. 
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Figure  A. 8:  Monthly  time-series  from  May,  1983. 
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Figure  A. 9:  Monthly  time-series  from  June,  1983. 
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ure  A. 10:  Monthly  time-series  from  July,  1983. 
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SUBROUTINE  DAYSCL{QI.QL.PQ,TAU,DSaTSC,USaiCON) 

THIS  SUBROUTINE  COMPUTES  THE  AMPLITUDE  OF  THE  DIURNAL 
CYCLE  FROM  THE  FOLLOWING  INPUT  DATA: 

Ql.  SOLAR  INSOLATION  MAXIMUM.  WATTS  PER  METER  SQUARED 
QL.  HEAT  LOSS.  WATTS  PER  METER  SQUARED  (USUALLY  NEGATIVE) 

PQ.  HALF  THE  TIME  INTERVAL  DURING  WHICH  Ql  +  QL  >  0.  SECONDS 
TAU.  WIND  STRESS.  PASCALS 

OUTPUT  VARIABLES  ARE 
DSC,  MINIMUM  TRAPPING  DEPTH,  METERS 

TSC,  DIURNAL  RANGE  OF  SEA  SURFACE  TEMPERATURE.  CENTIGRADE 

use.  AMPLITUDE  OF  THE  DIURNAL  JET,  METERS  PER  SECOND 

ICON.  FLAG  -  1  IF  CONVECTION  LIMIT  WAS  REACHED  (VANISHING  TAU) 

REFERENCE  TO  THEORY  IS  IN  PRICE,  WELLER  AND  PINKEL,  JGR,  1986 

DOCUMENTED  BY  J.  F.  PRICE.  JULY  2.  1985.  W.KO.I. 

COMMON/APLANE/BETA 1  .BETA2.aF>kLPHA,aCPW.RO 

THE  COMMON  APLANE  DELIVERS  PARAMETERS  WHICH  ARE  SITE  SPECIFIC: 

BETAI,  EXTINCTION  COEFFICIENT  FOR  LONGWAVE  INSOLATION.  METERS 
BEAT2.  AS  ABOVE  FOR  SHORTWAVE 

a  THE  FRACTION  OF  INSOLATION  ASCRIBED  TO  BETA  1,  NON-D 
F,  CORIOLIS  PARAMETEa  INVERSE  SECONDS 

ALPHA,  THERMAL  EXPANSION  COEFFICIENT,  KILOGRAMS  PER  METER  CUBED 
PER  DEGREE  CENTIGRADE  (IN  A  LINEAR  STATE  EQUATION) 
a  ACCELERATION  DUE  TO  GRAVITY,  METERS  PER  SECOND  SQUARED 
CPW,  HEAT  CAPACITY  OF  SEA  WATEa  JOULES  PER  KILOGRAMS  SECOND 
RO.  DENSITY  OF  SEA  WATEa  KILOGRAMS  PER  METER  CUBED  (CONSTANT) 


ICON  -  0 

COMPUTE  PS.  THE  ACCELERATION  TIME  SCALE 
PS  -  (SQRT(2.)/F)*SQRT(1.  -  COS(F»PQ)) 

QNET  -  Ql  +  QL 

CHECK  THAT  QNET  AND  PQ  ARE  »  0.  IF  NOT.  RETURN 

IF(QNET.LT.0-OaPQ.LT.0.1)  THEN 
TSC  -  0. 

use  -  0. 

DSC  -  999. 

ICON  -  9 
RETURN 
EMD  IF 


mmm 


CONST  -  (17RO)*SQRT{-ALPHA»G/CPW) 

C 

C  EVALUATE  TAB  SCALES  IN  THE  STRESS-DOMINATED  REGIME 
C 

use  -  15*SQRT{QNET*PQ)*CONST 

DSC  -  0.45*(17RO)*TAU*PS/(SQRT(QNET*PQ)*CONST) 


TSC  - 1 .5*(QNET*PQ)**  1 .5*CONST/{TAU*PS*CPW) 

C 

C  TAKE  ACCOUm-  OF  THE  EFFECT  OF  PENETRATING  INSOLATION 
C 

RSI  -  (1.  -  R)*(QI  -  QL)/QI 
HLAM  -  (1.  -  RS1‘EXP(-DSC/BETA2)) 

C 

TSC  -  TSC*HLAM«1.5 
use  -  USC*HLAM**0.5 
DSC  -  DSC/HLAM**1.5 
C 

C  NOW.  CHECK  TO  SEE  IF  CONVECTION  LIMIT  IS  REACHED 
C 

CALL  CDEP(Q(.QL.flBeTA  t.BETA2.CDZ.QOC,RfC) 

C 

TCON  -  PQ*QDC/(RO*CPW) 

C 

IF{TCOMLT.TSC)  THEN 
C 

C  IF  CONVECTION  LIMIT  WAS  REACHED.  THEN  USE  CONVECTION  SCALES 
C 

ICON  -  1 

TSC  -  TCON 

use  -  TAU*PS/(RO*CDZ) 

DSC  -  CDZ  +  (PQ/{R0*CPW))*R1C/TSC 
C 
C 

END  IF 
C 

RETURN 

END 

C 

C 

C 
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SUBROUTINE  CDEP(QI.QL.aB  1  ,B2,CD,QDC,RIC) 


THIS  SUBROUTINE  COMPUTES  THE  CONVECTION  DEPTH  FOR  THE 
DIURNAL  CYCLE 


IMMT  DATA  ARE 

Ql.  THE  MAXIMUM  SOLAR  INSOLATION,  WATTS  PER  METER  SQUARED 
QL.  HEAT  LOSS,  WATTS  PER  METER  SQUARED 
a  FRACTION  OF  INSOLATION  IN  LONG  WAVE  COMPONENT,  NON-D 
B1,  EXTINCTION  SCALE  FOR  LONG  WAVE  INSOLATION,  METERS 
B2,  EXTINCTION  SCALE  FOR  SHORT  WAVE  INSOLATION  METERS 


THE  OUTPUT  DATA  ARE 

CD,  THE  DAILY  MINIMUM  CONVECTION  DEPTH,  METERS 

QDC,  HEAT  FLUX  ABSORBED  ABOVE  CD,  WATTS  PER  METER  CUBED 

RIC,  THE  SCALAR  INSOLATION  AT  DEPTH  CD,  WATTS  PER  METER  SQUARED 


JIM  PRICE,  1  JULY  1985,  WKO.I. 


DZ  -  0.05 


SET  DEFAULT  VALUES 
CD  -  26. 

QDC  -  0.1 


DO  4  J-1,500 
Z  -  FLOAT(J)*DZ 

RIZ  -  Ql*(1.  -  (R*EXP(-Z/B1)  +  {1.-R)*EXP(-Z/B2))) 
DIDZ  -  QI*(R*EXP(-Z/B1)/B1  +  (1.-R)*EXP(-Z/B2)/B2) 
ELS  -  (raz  +  QL)/Z 


IF{ELS.GEDIDZ)  GO  TO  5 


CONTINUE 
GO  TO  9 
CONTINUE 


CD  -  Z 

QDC  -  (ELS  +  DIDZ)/2. 
RIC  -  RIZ 


CONTINUE 

RETURN 

END 
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